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ABSTRACT 


Compression testing of modem composite materials is affected 
by the manner in which the compressive load is introduced. Two 
such effects are Investigated in this report: (a) the constrained 

edge effect which prevents transverse expansion and is common to 
all compression testing in which the specimen is gripped in the 
fixture; and (b) non-uniform gripping which induces bending into the 
specimen. This study has developed an analytical model capable of 
quantifying these foregoing effects. The model is based upon the 
principle of minimum complementary energy. For pure compression, 
the stresses are approximated by Fourier series. For pure bending, 
the stresses are approximated by Legendre polynomials. 
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NOMENCLATURE 


<31. e 

n» £ 


®y 

Yxj^ or y ^2 


<1^. Oy or 03 ^. 02 
"y 


^xy» ^12 


■^xy* “^12 


xy* yx* 12* 21 


Eigenvalues 

Eigenvectors 

Strain In the x,y or 1,2 directions 

Shear strain In the x-y or 1-2 
plane 

Stress In the x,y or 1,2 directions 

Stress averaged across laminate 
thickness 

Shear stress In the x-y or 1-2 
plane 

Averaged shear stress 
Introduced material constant 

(ro"Ey/Gxy-2Vxy*Ey/Sx) 

Poisson’s ratio 

Poisson's ratio In the x-y, 1-2 
plane and their counterparts 

Angle between principal material 
coordinates and arbitrary body 
coordinates 

Ridgld body rotation 

Partial differential 

Complementary Energy Funoticr.nl 

C-ecmecric ratio ( Length, ■’’■■'id tn ' 


V 


det 



Ey or 





M, C, K or 
CM]I CC], CK] 

Q-I 

C-P 

U-D 

u« 


D«temilnan» 

7oung*9 acdxilus 

Young's aodulus In the x, y or 
I 9 2 directions 

Shear modulus 

Shear modulus In the x«y or 1-2* 
plane and their counterparts 

Matrices 

Quasi-isotropic 
Cross-ply , 

Unidirectional 

Rigid axial displacement 
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Chapter I 
33TERODUCTION 

A. Rationale 


In compression testing, it is difficult to 
determine purely compressive mechanical properties of 
fiber-reinforced matrix laminates. Some of the 
experimental data showed that mechanical properties of 
bhe specimen depend strongly upon the compression 
fixture utilised [!]• Therefore,- it is not surprising 
that some controversy has developed regarding acceptable 
techniques for compression testing. 

Other than manufacturing non-uniformities in 
test specimens, compression data may be suspect due to 
uneven gripping of the tabs, poor alignment of the test 
machine and/or poor alignment of the test fixture. 
Fracture or ultimate compressive stress may be difficult 
to obtain because another mode of failure (l.e. buckli.ng-, 
delamir.ation) may occur first. 

Many of the foregoing difficulties can be lessen- 
ed by selecting a short gage length specimen. At the 
extrenlties cf the test specire’^. the clamti.nc r.ec h.ar.isr. 



of the test fixture creates a complicated stress state 
by preventing transverse expansion. For sufficiently 
short gage lengths this constrained edge effect will be 
evident throughout the entire specimen. A size change 
of the specimen, therefore, may merely substitute one 
difficulty for another.' 

The constrained edge effect has been pointed 
out in Refs. [2-^3. However, only the present work 
provides a model capable of quantifying it. This is 
done through stress analysis by assuming perfect align- 
ment and two different gripping mechanisms of the 
fixture: (a) uniform gripping (axial compression) and 
(b) small in-plane bending superimposed upon axial 
compression. 


B. Background 


Pagano and Kalpln [2] investigated the Influenc 
of the end constraint, both experimentally and 
analytically, in tension tests of anisotropic bodies, 
including on off-angle graphite /epoxy laminate. Their 
analysis v.*as based upon the two-dimensional elastic 
conpaclbillry equations. They concluded that the gripe-' 


r.eohanisr. and the length to width ratio of the speci.ne: 
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were the principal reasons for the non-existence of a 
uniform stress state. Ihey also predicted a more serious 
influence of gripping in compression and torsion testing 
of anisotropic bodies. However, they did not quantify 
the end constraint. A photoelastie study of axially 
compressed rectangular sections, by Phillips and Mantel 
[3l> gave some evidence of the effect of load intro- 
duction upon homogeneous, isotropic materials. 

An investigation of the effect of an end attach- 
ment on the strength of fiber-reinforced axisymmetric 
coaqposite cylinders was presented by Whitney, Grimes and 
Prancis . They pointed out that an end attachment 
which allows some deformation of the end (e.g. adhesive 
bond) will help alleviate the problem of high stress 
and strain concentration at the attachment end. 

Another method of studying the edge effect in 
two dimensional stress analysis is based upon the Airy 
function. Unfortunately, the mixed form of the boundary 
conditions precludes any possibility of an exact solution. 
Hess [53 uses separable forms of the Airy function (which 
decay exponentially from the fixed end) to determine 
approximate solutions. 

A related problem, v.'hose solution also 



mor« attention In recent years, is the free edge effect. 
At the traction-free edge of a compression specimen, the 
mismatch in the material properties at laminate inter- 
face causes a hl^ly localized effect. 

An ezaa^le of the trtB edge effect fo:^ a biaxial 
stress state using methods developed in the present study 
would seem an interesting challenge. 

The difficulty in estimating stresses in the area 
near the free edge, using the finite difference method 
presented by Pipes and Pagano C ^ 3 was pointed out by 
Wang and Diclcson £9 ]. The finite element procedure 
developed by Wang and Crossman £7 ] has the same 
difficulties as the finite difference method. Both 
methods ne^ certain artificial manipulations, specifi- 
cally in the region very close to the free edge. The 
perturbation technique applied by Hsu and Herakovich £ ^ ] 
provided smooth continuous stress distributions in the 
vicinity of the free edge and mathematical evidence of 
singular ..iter laminar shear stresses for cross-ply 
graphite/ epoxy laminates. Another method of estimating 
che lr.cerla.T.lnar shear stresses is based upon the 
lalerhln method [ 9]. 

In the subsequent chapters, an analytical 
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:• Ir.g the ccnscralned edge 
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effect. 

I 

C. Specific Objectives 

The objectives of this study are: 

1. To develop a closed fora approzimation to the stress 
distribution within each lamina of high-strenf^h 
graphite/epoxy during compression tests. 

2. To determine the effect of specimen geometry upon the 
measured compressive properties, including the 
determination of the minimum specimen gage length 
necessary for the existence of a uniform compressive 
state in the central region of the specimen. 

3. To determine the effect of small in->plane bending 
upon the meastnred compressive properties. 

4. To determine the optimal location of strain gages 
for compression tests. 


Chapter XI 
PROBLEM FORMULATION 

A« Statement of the Problem 

The primary objective of this study Is to 
determine the effect of testing devices on the response 
of compression (tensile) specimens of laminate composites 
which are symmetric about their middle plane. For the 
case of perfect alignment and perfect gripping In a rigid 
flxtiire, the ends of the specimen will undergo the rigid 
displacement shown In Pig. 1, where u* and 6® denote the 
unlfoxm displacement and rotation of the constrained 
edges » respectively. Also, L and b are the respective 
half length and half width of the specimen, and x and y 
are Cartesian coordinates measured from the specimen's 
center. 


B. Kodelling Assumptions 


syrmierrl 




The laminate theory for fiber layups which 
c about the riddle plane Is applied here, 
u" zhe laminate ohere exists a generalized 


are 

Thus 
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state of stress whose Cartesian components are denoted by 

bars above the stress symbols Indicate 
quantities averaged across the laminate thickness. 

With the assumptions of small displacements and 
a linear orthotropic constitutive response, the field 
tiquatlons to he satisfied are: 

1. ecrulllbrlum equations 

+ 3x/3y - 0, 

8t/3x + &ff_/ay - d; (1) 

, . y 1 

2. strain-displacement relations 


( 2 ) 


( 3 ) 


■ au/3x, 

Sy - 3v/3y, 

Y • 3u/3y + 3v/3x; 

3. constitutive equations 


G » 

S . , 5 + 


X. 

11 X 

12 y* 




■ 


®22‘'y> 

y “ 

^ 44 ' ; 
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In Eqs. (2)» u and v denote the displacements in the x 

« 

(loading) and y (transverse in the plane of. the specimen) 
directions, rdspeotively. The material constants Sij can 
be computed directly from the known material constants of 
the constituent laminae and their fiber orientations with 
respect to the x-axis p.03. 

The boundary conditions which are to be adjoined 
to Eqs, (1-3) are of mixed type. On the stress-free edges 
we have the static boundary condition 

Oy • T - 0, on y - + b. 

On the other hand, the kinematic boundary 
conditions, according to Pig. 1, are 


u(+L,y) "^(u® + 0®y). C^'j) 

D\ie to the linearity in the Eqs. (1-3) and 
boundary conditions (4b), the Principle of Superposition 
is applicable and it suffices to solve the purely 
compressive case (9®»0) and the pure In-plane bending 
(u®«0) case, separately. 

A common method of obtaining the solution to Eqs. 
(1-3) is based upon the Airy stress function Cll]. The 
resulting fourth order equation is generally solved by 
separation of variables. The mixed boundary conditions 


(la, b) make this approach extremely cumberscm.e. 
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An alternative approach Involves reformulating 
the problem In terms of the complementary energy [7,8] 

. ♦ • f CSu'i “l2Vy * 

+S|jljT^]dx*dy 

“ *^-b Cu(L,y)a^(L,y) 

+u(-L,y)a^(-L,y)]dy. (5) 

The Principle of Minimum Complementary Energy 

states that of all stress fields (o_, ? » t) that 

X y 

satisfy the equilibrium equations (1) and the static 
boundary conditions (4a), the exact solution actually 
minimize Thus Eqs. (1<>4) may be expressed simply 

' I 

64 • 0. X 



Plg.l Edge displacements due to perfect 
gripping and perfect alignment. 


Chapter III 
FOBHAL SOLUTION 


A. Pure Compression 

Since, in this case, 8°-0, It Is evident that 
axial stress must be an even function of y, and that 
the shear force along the edges, x*±b, must vanish. It Is 
therefore clear that the stress may be represented by 
a Fourier series. 

cos^ 

For N sufficiently large, and for fixed x, the series 
will uniformly approximate on all .Intervals for which 
dj^ ls continuous. At points of discontinuity for the 
series converges to the average value of a . Here, 

(n»l,...N) are the unknown Fourier coefficients, 
and c Is the average compressive stress across every 

w 

x»ccnstant Section. 

The remaining stresses a . and x are obtained by 

y 

solving the equilibrium equations (1), using the boundary 
condition (^a). These results are 


®7 ■ ”4 * (- 1 )”*^} ( 6 b) 

r-,^aln!?Ip;(If). (6e) 

In deriving Sq* C^c)« the vanishing of shear 
force- V •; 

- 0 . 


was used to detemlne the constant of Integration* Also 
(, y indicates differentiation with respect to indicated 
argument. 

The unknown functions Pn niay be determined by 
substituting Eqs.(6) into Eq. (5)» that is 




-Cl+ r n ooaSgi P„C^)] 
+2Sj^ 2C1+ £ m oos2^ 


M . 

+ E - [cos 


rain’’ 


'22 
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N 
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1 r nuy 

- Lcos-^ 
n D 


+ (- 





n'' b' 
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K N 

+S,4C I r slnSsZiF^(I|)]j(lX(l7 

* 1 X 

u 

Cl+ t n cosSIi P (irC)] 
e ~u 2^ on 


N 

+Cl+ E n cos 
1 




(T) 


whe^e E*L/b . 

After Intergratlng over the y coordinate, appli- 
cation of standard techniques of Variational Calculus 
[ 13 ] renders the expression 


Si 




n 


+S22 


2^ n 0 • tt*n m 


=^ 1.4 f 

12 ^ ^ f —1 ^ 

+ K ^-^12 I '’n<±^5)«Si2 2 


vn 


+3S22I <(i'C )*2S22^J=^, 


•p'(+»c)]-cs?l(+ic)] 

Q • ** *^ 


15 N , N , II, 

+ [S,, I ?^(+’:a- 3S,,i: ^ ?^(+ttO 

<: c 12 . n £2, „2 n - 
1 in 


-2s 


izil 


m+n 


w 


- — P (+trC ) 

22n»in “ “ 

Nr 

+ Pj^(+n)]'C«P^(fr01 . 

- 0 . ( 8 ) 

Xt Immediately follows from (8) that the Eixler 
equations are 


N 


I .CS222jpW(I|). (2 s^ 2«^4) P;(If)+Sj2 P„(2|)] 


>m+n 


* ^ Pf (If)] 


• 0 . 


and the natural boundary conditions are 


(9a) 


X n 


m4n 


+ P„(i^)] 

« r [2S22(-l)”'^Vn], 

N 

I C(S„+Sj^2) P^t+’fO-Sjs^ Fn(i»C>: 
1 - • 

/ T ttr 

-r ■■ F (~-r)] 

22 n*n r. - 


(9b) 


ni3»n 


»0. 


(9c) 


Ml- 


With the introduction of the followint; definitions 


S<« ^ ■ 1/E • S, - ■ -V /E • -V /E » 
11 X* 12 xy X yx X* 


Ijj i - I/O 

r, - C2S„ + s^^j/s 


K • 


12 “M"“22 

3 2,-1, 2,+l, . 

7z - * 


xy* 

E_/a_ - 2v 


i 

22 


Sym. 


y' ’xy - yx* 


2f-l^ + 

jC-j) + 

a . . 

32 


1 

n2 



E /E 

y X 


9 



Eqs. (9a, b, e) may be conveniently cast in the 

form 

M - To 1 •»- K ? - 0, (10a) 

mm •• 

and 

5 ^ ■*• V - 

Cfo + Vj^) 1 - M P*(+vO - 0. (10c). 

It is also necessary to determine the constant 

After differentiating • (Eq, 7) with respect to o^ 

• • 

and simplifying the resulting expression with the aid of 
Sqs.(9a,b,c) and setting 11 ■ o 

* 

c 

we obtain 

N , - t 

^ - # s 5 

It should be noted that if Pxi*0 (n*l, ...N), then Eq. (11) 
reduces to the elementary strength of materials formula 
[ 12 ]. 

The solutions for N; terms retained in the series 
(6), and the solutions to Eqs. (9a,. b, c) give the "best" 
(In the mean square sense) N-tj'm approximation to the 
true solution. Therefore, It is reasonable to expect 
that the approximate stresses vrill be closest to their 
true values at locations v;here or.elr 


true \'alues are the 


largest, that is, at the constrained edges. 

Eq» (ina) is a fourth order ordinary differential 
equation with ‘constant coefficients. Thus, we assune a 

solution or the fona 

• • 

F ■ XL cosh(oir3c/b). 

Substitution of Eq* (12) into (10a) generates the 
eigenvalue ^ prob lem 


Co*M-o*r.l + K] n - 0 


(13) 


for the eigenvalues n and eigenvectors n. 

Note that a and ^ are obta^ed independent, of the 
boundary conditions, and hence they do not depend upon 
the manner in which the compression load is introduced. 
They depend. only upon the number of terms N retained in 
the series, and the material constant Fq. in the case of 
quasl-isotropio lajrups, Fo«2, and they are 

also Independent of the material constants. 

A necessary and sufficient condition for the 
existence of a non-trlvlal solution to Eq. (13) is 

detCa^ K - a^Fo 1 + H • 0 Cl4) 


Eq. (14) is a polynomial of order ^N; however all 
solutions must occur in equal and opposite pairs. And if 
the roots are complex, they must also occur in pairs of 
complex conjugates. 
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If, for complex eigenvalues, the components of 
the p-th eigenvector are assumed as 
2p Ap Cl Up,- • • Hp„]. 

then the general solution to £q. (10a) for even 

functions becomes 
n 

N 

P- ■ I {A„ n „cosh(o^vx/b) 

" p-1 P P” P 

+3(_ Tj__cosh(o vx/b) 
p pn p 


N 

» 2«Re{ Z k n coahCa vx/b) (15a) 

p.l P Tn p • 

where bars above the symbols indicate complex conjugate. 

1 2 

For real eigenvalues a*, o^ (p«l,2,* • • N) 
the general solution Is 


N 1 1 1 

P ■ E [A" cosh(o“vx/b) 
n P pn p 

+ Ap cosh(apTTx/b)3, (15b) 


where 


3 •"'* 


«1 - 41 
-pn ^p 

Cl 



-a 2 
-pn ? 

Cl 

^p2 


The ccr.pl 

ex 

constants 

A or real 

P 

■e decerri 

ned 

, fror the 

tcundary c 


D 


lot, C) 
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The remaining possibility that some roots are 
i*eal and some complex was. not encountered In the 
numerical calculations for the assumed data* Thus* 
although this csise Is.^ as routine as the two foregoing 
eases, it will not be discussed further. 


B. Pure in-plane Bending 


For pure b'sndlng* uo«0* and the axial stress 
is an odd fimctior of y. As before* the shear force 
along 'the edges x*±b must vanish. However* it turns out 
that a Fourier Sine: series approximation to is not 
convenient. This is because application of the boxindary 
conditions (4a> to the stresses obtained after integrat- 
ing Eqs. CD introduces side constraints on the Fourier 
coefficients. ' In order to circumvent this difficulty* 
is expanded in terms of odd Legendre polynomials. 

One definition of the N-th Legendre polynomial 


la C13] 


V C2n-2k)! r 7 >n-2k 

fc-0 2^1 Cn-k)! Cn-2k)I ^ 


where N * 


n: even 




Cn-1) 


n : odd . 


C 


(16) 
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Th«n the stress cay be represented by . 

'.-“b Vf>- 

Here, u*e the unknown generalized Fourier 

eoefflelenta^ and la a constant. . . 

The remaining stresses o . and t are obtained by 

JT 

solving the equilibrium equations (1)* subject to the • 
boundary condition* (4a) . These results are 

* - - ? r 

"y - (2n-lKJnf3) 




(17b) 


N 


Wb’ ■ 

m deriving Eq.' Cl7c)« the vanishing of shear 


force V 

•^-b - Q 

was used to determine the constant of Integration. Also 

I 

( ) denotes differentiation with respect to the indi- 
cated argucent . 

The unknown functions may be determined by 

n 

substituting L-a. (17) into Sq. (5), and Integrating 
over the y-cocisilnate. After applying the variational 
nethoc and collecting terns, the Euler equatlcns are 


- c t alne c • 


( 18 ) 


whsr« 


M - To C Q + K 

m *- m m 


a - 0 


MCN.1I) 2 , 

(2M-3) (2»-l) {2m) (2N+5) 


M(N-1,N) - 


(2N-5) (a^-3) (2n-l) (2N+1) (2N+3) 


m ( m ; n ) 


(2N-7) (2If-5) (2N-3) (2M-1) (2Mfl) 




2 . -1 
9-11-13 11*13*1S 


L\^l) C2I'I-i 



•» 



K - 
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^ 0 

• 

■V=x -2- 

°3 

°5 
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'e 


* 0 
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(2ii+l) 




The natizral boundairy conditions associated with 
£q. (18) are 

m g" (±5) + C G (±C) - ''yxC 10 0. . 0]^ , (19a) 

xy - - 

M g'" (±e) - (ro+v,,)C Q* (±5) » 0 . (19b) 

- — . yx - — • 


d4 

In addition, the condition 5J”*0 results In, 

b 

after considerable manipulations. 


to . E 


jt (03 . 


E^.2b 


( 20 ) 


where 


“b “■•'■-b y ®x * “T" "b 


as In 


It should be noted that In Bernoulll-Suler theory 

8® - - A - (21) 

-x*^^ 

The method of solving Eos. (Ir, 19) is the sane 
:he pure compression prctlem. write 


a-icosh(^) , 


where g and £ are the eigenvalues and eigenvectors , 
respeotivelj* * Thus 0 and ^ are determined from 


M - 0^ FoC -I- KI i - 0 . 


( 22 ) 


If the components of the complex eigenvectors are 

Sp-S&S2 • • • 

the solution for 6n ia given by 

For real eigenvalues « the counterpart to Eq. (15b) is 




(23b) 


C. Lamina Stress 


In the previous two sections (A, B), we formally 
obtained approximate solutions for the average '.stress arsd as. 
a. consequence of Eqs* (3)» strain. The remaining: task Is 
now to obtain the stresses within each constituent lamina. 
The approach here will follow Jones Cl03* 

For either quasi-lsotropic or cross-ply 
laminates, the stress-strain relations (3) nay be written 




^ 1 -V 

f E" 0 


m • 

“* 


^ m 



► - . 

r 1 ° 
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• 'y 

' - CQ3- 

“y 

' (24) 
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to 
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wher« E la the Young's modulud (E^* Ey* E) , v la the 

Poisson's ratio (v. ■ v • v), and Q la the shear modulus 

(0 • G * G ). Each constituent lamina of the laminate 
sustalpa the same strain Ce^, Sy» 3^ m the x-y 
coordinate system. 

However j the stresses differ from one lamina to 
the next. It Is necessary, therefore, to determine the' 
appropriate stress-straln relation for each lamina.' Let 
us suppose t^t the principal material axes are Inclined 
at an an^e e to the x-axls (see Fig. 2 ). Then the 
strains In the material coordinates are obtained from 
the laminate strain by 


- . 




’*2 2 

^1 

k." CT] , 
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C(2] - 

cos 6 sin 6 sln26 
2 2 

sin 9 cos 6 -sin29 

^12 
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. 2j 


-sin20 sin20 cos2» 
2 2 


( 25 ) 


Similarly, the stresses In the principal coordinates are 
given by 


C261 


Row let the constitutive equation in the 1*2 
principal coordinate system be 
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^12 
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(27) 
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where Ex, E2 are lounge’s moduli in the 1 and 2 directions j 
respectivelyjVj^ia Polsson^s ratio for stresses applied 
in the fiber direction ^ 
principal shear modulus. 

iUrter combining Eqs* (2^1-, 25, 26, 27), we obtain 




m ^ 

0 

X 



“y 

► - [Tr^CS 3 [T][Q] ‘ 
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X 




( 2 °) 


where 



r ■ 2 

cos 9 

sln^9 

•sln29 

Ct3“^ - 

sln^e 

2 

cos 9 

sin2e 


sinS 

~sin6 

cos29 






Following fundamental matrix algebra, Eq. (28) 
can be simplified to the form 




a 
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V 

1—1 

0 

i_j 

1 

A 

"y 



T 

V y 


where 

Ki - CTr^-cs:-m-[Q]. 

CXeairly,. the matrix [C] depends upon the lamina 
material properties and orientation* 



Chapter IV 






NUMERICAL SOLUTION TO TEE EIQEK7ALUE PROBLEM 


We now take up the numerical solution to the 
polynomial of order 4 n 

f(a) • det Co^M - o^FoC + K] - 0 • 

Clearly f(a) has the form 

f(o) + C 2 ,a^^“^ + ... + Cjjo^ + ^ 

where Cj^ are functions of certain invariants of the 
matrices M, C and K and their products. For example, 
Cjj+i^et CKl and Ci*det Dl] . Ths other coefficients, 
however, are considerably more Involved. 

A numerical method for determining the set of 
Cl’s for given N relies on the utilization of a high 
capacity computer. By reasonably choosing a. set of 

t T » 

arbitrary numbers 02 » • • • ®n+i 5 evaluating 

f ■ ■ * 

f( ) N-<*1 times, we obtain the simultaneous equations 

f(aO 


C, (a')'*" * C2 

* ^2 f“2> 


. . + c 


N+1 




+ . 


T hlj f htf—O t 

The ooefflelents Cji, Basr now be- routinely 

obtained by solving the above set of simultaneous 

equations In which C^'s are the unknowns. The numerical 

sensitivity of the procedure may be checked by choosing 

several different sets {a^> and comparing the solutions. 

Furthermore, in our ease, it is convenient to reduce the 

» 1/2 

order of the polynomial from 4 n to 2K, by taking (a^ )■ 

I ’ 

Instead of . This step also speeds up the process 

I 

of obtaining the roots of the polynomial. These roots, 
i.fr.. the eigenvalues, were obtained by using a standard 
sxzbroutine based on the Newton-Raphson method. [14]. 

^genvecturs . . . are determined from a set of 
linear algebraic equations ( 13 or 22 )» Since the 
eigenvectors are not unique, a very convenient 
normalization procedure is to set the first component of 
each eigenvector equal to unity. The coefficients Ap, 
determined by the natural boundary conditions, are also 
routinely obtained by solving a sot of linear algebraic 
equations . 


• CWiptei* 7 
NUMEBICAL RESULTS 

A> Pure Conpreasion 

■ a. Quasl-Isotroplc (Q-I) C0/±45/90]3 Laminates 
X- Generalized Plane Stress 

t 

Aeearding to elementary rod theory « the stresses « 
sufficiently far from the edges at which the load is 
instroduced, are assumed to be uniaxial, i.e. 



Thus the specimen must have a sufficiently long gage 
length if Eq» ( ) is to be applicable anywhere. 
Elementary teecry/- however, is able to provide neither 
the minimiun gage length necessary for Eq. ( 30 ) to hold 
nor the stresses in the neighborhood of the clamped 
edges . 

In the present approach, the general plane 
stresses ax> ^ «i®P®nd upon the material constants 
Vyj. , Ey/Ej^ and To, and the specimen geometry ratio C. 




For quasl-lsotropic laminates, ro»2, » and hence 

A y 

the stress distributions merely depend upon v and C* 

The stresses were calculated for a range of C 
from 1/k to 6, and for v«0.336 which is a fairly typical, 
value for graphite/epoxy quasl-lsotropic laminates. It 
will be noted from Eqs.(lO) that the stresses are 
approximately propotlonal to v and therefore approxi- 
mate solution for other Poisson ratio *^s may be obtained 
by sealing the ctirrent "solution. 

Since re>2 for all Q-I laminates, the eigenvalues 
(1"1,2, • • • ,N) computed from Eq.(l4) merely depend upon 
the number of terms, N, retained in the series [Eos. (6)]. 
The results of this computation, as explsined in Ch. HI, 
are shown in Table 1 for values of N ranging from 1 to 
10. It will be noted that all the eigenvalues for.n<i0 
are complex values, and consequently the solutions for 
the functions (n«l,2, * * • ,N) are given by Eq. (15). 

The stresses [Eqs. (6)] were plotted for 
different values of N at various cross-sections x/L ■ 
constant in order to assess convergence for Increasing N. 

* This values was obtained for the material properties: 
Et“21x 10^ ksi, E2*1-7 x 10^ ksi, 

G^2®C*^5x 10^ ksi, E*7. 95x10^ ksi. 




0.59 0.67 

9.65 
































































































Aa shown In Figs . (.4,6-8) , convergence was excellent 
•way from the constrained edges with just three terms 
retained in the series.. Raarv'bat: not at the* e lamped 
edges; coxivergenee' was excellent with, only six terms^^ ' 
retained In the FouTier Tserles 'solutions Esee Fig. 5v 
14].. Another measure of the - cotrrergence is' provided by 
Eq. (11) . In Table .2 the -ratio: was evaluated 

for various values of N and 

The convergence of for Increasing N is 

evident from Table 2. Note that f-sr large the value 
approaches unity, Milch is the result predicted from 
elementary rod theory, Eq. (30). ©le reciprocal of the 
entries in Table 2 represents the apparent percentage 
increase in average stiffness due to the constrained 
edges . 

It is convenient to write the generalized plane 
stresses in the form 

?3j*-CgCl+6jjCx,y)], 

Oy* o-^,6y(x,y), (31) 

Clearly, for |6y| and (6^y| sufficiently 

small, Eq. (31) ..'-'ll closely aoproxinate Sq. (30). We 
shall say that the stress state Is aoorc.’clr.ateiy uniaxla 


TABLE 2. values of u^E^/o^L for Q-I laminate 
In pure compression 


. ;u«Ej/a^L 


0.25 0.50 


0.9342 


0.9203 


0.9173 


\ 

3 

.6 

0.9926 

0.9963 

0.9907 

0.9954 

0.9903 

0.9952 

0.9902 

0.9951 

0.9902 

0.9951 




























33 


at a given x*conatant cross-section if 

I < 0.02 

I ;^0-02 for |7 I ^ b ^32) 

< 0.02 

The 2 % bound on the deviation of the true stresses froa 
the uniaxial state Is, although arbitrary, quite useful 
particularly for the experlaentalist . By providing a 
definite bound, the effect of the constrained edge can 
be quantified. 

For geometry ratio's C<1*5, (32) was not satisfied 
anywhere. Thus the effect of the constrained edges is 
observed everywhere in the speclaen. Stresses distribu- 
tions along the center line x**0 and at the edge x*L are 
shown in Pigs. 6^-11 for C*l/^> C“l/2, and C*l» respectivelv. 

For C^l«5, there exists a region in which the 
stress state is approximately uniaxial, l.e. Eq. (32) 
is satisfied. It was found that the domain of influence 
of the edge is limited to 1.5b Cor 75S of the width), 
as depicted in Pig. 3. 

As expected, the stresses in the shaded region 
of Pig. 3 are Independent of C provided L>1.5b. 
Consequently, once the stresses are determined for one 


value of C> 1 . 5 . 


are 


>nr--.T. f:r all values 




Influenced by 






Increasing L, for fixed b, merely Increases the uniaxial 
stress domain. Results are shown in Figs. (4,5, 12> for 
C=»3 and x/L» 0;T5> 0-9 and l.p respectively. According 
to the foregoing discussion, the generallxed stresses 
are the same for C«6 and x/L»0.875, 0.95 and 1, respec- 
tive y CFlga. 13-15] 

It should be observed that the stresses at the 
edge x*L for [Fig. 11] and 5*3 [Fig. 12] are almost 
Identical. The reason for this is because the stresses 
at x»L are affected by the constrained condition at only 
that edge; the stresses at each edge x»±L for both 
and 5«3 are outside the domain of Influence of the other 
edge x»+L. 

As Indicated earlier, the stresses at the clamped 
edge appeared to be converging CFig.12] quite well for 
N*6. A closer examination of the tabulated values of 
Ox(L,y) did indeed confirm convergence for |y|<b.^ 
However, at the corners y»b, the stress appears to 
grow without bound. This apparent singularity is shown 
in Table 3- 


2. Lanina Stress 


Figures l6-i8 show representative laminae stress 


at z'r.e conarrained edge of C-I laminate 


^ ^ 


c c L i *8 s 






















































STRESS DISTRIDUTIOM DIAaRAN 



PlKure l| Compression, Q-I laminate, C-3»x/L»0.75, M*3,H 










PlRure 5 ConipreBBlon, Q-I laminate 






STHESS DISTRIBUTIOH DIAQRAM 


I 








/b y/b y/b 

Figure 6 Compreaslon, O.-l laminate, x/L“0, N»3»*l 








STRESS DISTRIBUTION DIAORAH 



Figure 8 Compreaslon, O-I laminate, C"l» 3t/I«*0t N*3**l 








STRESS DISTRXBUTIOM. DIAQRAM 



PlKUre 10 Compreaelon, C'-I laminate, C*l/2, M»6,7 







STBESS DISTRIBUTION DIAORAM 
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ic 


Figure 11 Compres 8 ion» Q-I laminate, C* I » x/L»l, 









STRBSa DISTRIBUTION DIAORAM 



CompreBslona Q~I laminate, t*6»x/L"0.fl75* 











flTRBSa 



ComproBBlon^ O-I laminate, C-6, x/L-0,95, N«*6,7 








STRBSa DISTRIBUTIOM DIAQRAM 
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15 CompresBlonf Q-I lAmlnate, N*«6,7 




STRESS DISTRIBUTION DIAORAM 



Compression • Q-I lamina 





STRESS DISTRIBilTIOM DIAdRAM 









Figure 37 CompresBlon 







STRESS DISTRIBUTION DIAORAM 



Figure l8 Compression, Q-I lamina 







oriented at‘0-0®, 45®, 90® respectively. The particular 
plots are for C-1/4. 

b. Crcfss-Ply (C-^) CtUS]^ Laminates 

Oen«*alized Plane Stress 

Unlike Q-I laminates, C-P laminates have a nei^-> 
tlve material constant ro»-1.77*f and a relatively large 
value for Poisson's ratio (v«0.801 ). As a result of the 
hi^ value of v, the influence of the constrained edge 
should be expected to be much greater than for Q-I 
laminates, 

. Just.- like for Q-I laminates:, • the eigenvalues are 
again complex. [See Table 4]. Convergence was somewhat 
slower than for Q-I laminates;, more terms were needed to 
obtain a reasonable approximation to the stress- distri-’ • 
butlons at the edges. Plgiires 19 and cO show the stress 
distributions at x*L, for C*l/4 and E»3, respectively, 
for N-9 and 10. However, for the region x/L<0.9, 
a 6 -tera approxliiiatlon showed excellent convergence. For 
example, see Fig. 21 for which C"3; stresses at x/L»0.9 
are plotted for N=6 and N-T. 

* This value "was obtained forT 

E]_»21xlQ^ ksi , E2=1.Tx10- ksi , V]_2=0.21 , 

Gi 2*0 . 65x10- ksi , Z“T.9-jxl0^ ksi 



8.72 
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i 

, _ 1 

1 

Once again there appears to be a singularity at ! 

the comers y*±b. The data, tabulated in Table S, ! 

c^tainly do not suggest conrergence. 

i 

In Table 6, the values of u®E_/a*I« are 

• • . X c 

tabulated for C-P laisinates, and the conver^nce for 
increasing N is slower. Also, observe that for very 
small aspect ratios, the apparent stiffness increase 
is well over lOOS. 

According to the aiforedefined axial stress state 
CEq. (32)3, the domain Influenced by constrained edge 
turns out to be precisely double that of Q~I laminates. 

Consequently, an aspect ratio fbr which 5-3 is the 
smallest length-width ratio for which Eq. (32) is 
satisfied along the center line x*0. For C>'3 » a 'uilaxial 
stress field will exist in a region '-around the center 
line x«0; the range of length 3b, the domain influenced 
^7 the constraint has length 3b, measured from the edges 
x*+L. 

2. Lamina stress 

Since the lamina s-cresses are linear combinations 
of the averaged stresses o’ , -J , 7, [see Eq. (29) ]j lamina 
stresses will converge at the ssme rate as the average 
The stresses at the corners ?.c the clamoed 


stresses . 
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TABLE 5. Normalized stress at corners o^{L,^)/Cq 
for C-P laminate in pure compression 



10 


1.3^ 


2.09 


3.^2 


3.26 



















































TABLE 6. Values of uPEj^/a L for C-P laalnate 


In pure compression 




0.25 0.50 


3 

6 

0.9335 

0.9668 

0.9183 

0.9592 

0.9147 

0.9573 

« 

0.9137 

0.9569 

0.9131 

0.9565 





































STRESS DISTRIBUTION DIAORAM 



l'lEure?2 CompreBslon, C-P lamina 
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edges will exhibit sigulari ties. However, for the region 

away from the edges, the stress distributions appear 

% 

very well behaved. For example » see Pig. (22) in which 
x»0 (at the center line), N«10, 0»45®. 

c. Unidirectional (U-D) CO]. Laminates 

B 

Since all fibers lie in the same direction, the 
generalized plane stresses and lamina -tresses are the 
same. Also 2y"E2, and 

assumed data Q^*0.65xl0^ ksi, v^*0.21, we compute 
ro-2.58l. 

' Unlike the previous two laminates, the eigenvalues 
of U-D laminates were real. Table 7 lists the eigen- 
values for up to 3 terms. Figure (230 shows the resulting 
stress distributions at the edge for C*3» The values of 
stress 5^ at the corners for different C are tabulated 
In Table 8 for N ranging from 1 to 3* According to our 
definition, a uniaxial stress state does exist everywhere 
except at the corners. 

















STRESS DISTRIBUTION DIAQRAM 
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Compredslon, U*-D laminate, x/L"l, N»2,3 
























B. Pure Bending 


Conversence of the solution in bending was faster 
• * 

than for pure compression for both the Q>I and C-P 
laminates » This suggested that the Legendre polynomial 
may be preferable to a Fourier series for similar mlzed 
boundary value problems. Since the axial stress o^. when 
reduced to elementary bending stress. Is linear In y. 

It is apparent that Fourier series will take many terns 
to approximate cr^ In y-coordlnate while the first order 
Lege..dre polynomial Is equal to y. 

a. Quasi-Isotroplc 
1. Generalized Plane Stress 

Eigenvalues are shown in Table 9 for N ranging 
from 1 to 8. Observe that the eigenvalues are increasing 
at a faster rate with the number* of terms for N>5 than 
the corresponding case in pure compression [see Table 1]. 
This may account for the apparent faster convergence of 
these stresses. Fewer terms for approximation of the 
stresses are needed than for pure compression [compare 
a in Figs. 12 and 24]. 

The effect of the constrained edge is comparable 
to the pure compression case. Outside of the region of 
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laminate, 5 x/L«l, H*»5,6 







Influence of the constrained edge, the axial stress Is 
linear and the other stresses vanish [Flg.25] exactly as 
predicted by elementary theory. 

Another measure of the constrained- edge effect 
is provided by Table 10, In which values for 
0°2b^/3J%Sij^C are tabulated for various C ahd N. 

Again, note. that as N Increases, the value rapidly 

converges. Also as? increases, the value 0®2b 
tends toward unity, the predicted value from Bernoulll- 
Euler deflection theory. 

•A possible stress singularity at the corners, 
x*±L, y**±b. Is very much in evidence from the stress 
plots In Pig. 2U. Alternatively, the value of a-^(L,b) 
Is tabulated for different M in Table 11. and shows no 
sign, of converging.. 

2. Lamina Stress 

Figures 25-28 show the laminae stresses at the 
constrained edge for the plies oriented at 45®, 

90°, respectively. Xt is Interesting to observe from 
?ig. 28 that the greatest normal stress occurs in the 
direction of the fibers (i.e. the y-dlrection) . 

However, in the 0° lamina (Fig. 26 } the greatest stress, 
except for the corners, occurs transverse to the fibers 


























































STRESS DISTRIBUTION DIAORAM 



Figure 25 Rending « Q’"I laminate, x/L“0,5, N^SfC 
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lamina 
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(also the y-directlon) . 


b. Cross-Ply Laalnates 
1. Generalized Plane Stress 

The eigenvalues for pure bendlna were once again 
complex valued and are tabulated in Table 12. The 
accelerating rate of Increase of the elf^envalues is 
apparent from the Table, and is reflected in the rate of 
convergence of the stresses [see Pig. 29 j. 

Table 13 provides values for 6®2b*/3M^S3^3^C 
/ ^ 
and is the counterpart • to Table 10 for Q-I laminates. 

As we observed for pure compression, the effect of the 

constrained edge upon apparent bending stiffness is 

considerably greater for C-? laminate than Q-I laminate 

In pure bending also. 

Again, evidence of a singular stress state at 
the corners of the clamped edges is provided by Fig. 29. 
Tabluated values of o^^C+L, ±b), as shown in Table l^i., 
also appear to grow without bound for large N. 

2. Lamina Stress 
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TAPLE X2, EJgenvaineq fpr CrP laminate In pure bending 































































































































c. Unidirectional (U-D) C0]_ Laminates 

s 

Flj^ur'e 32 shows the axial stress Is linear and' 
the other stresses vanish,, for 5*6, x»+L. exactly as 
predicted by elementary theory. Eigenvalues are real and 
are tabulated In Table 15, for up to three terms. Again, 
calculations for axial stress at the corners of the 
clamped edge suggest a possible singularity. 






















Chapter VI 
SIGNIFICANCE 


A general method of solvlzig the two dimensional 
stress analysis problem for rectangular laminates 
subject to mixed boundary conditions has been presented. 
For compression sepclmens, the kinematic boundary 
conditions define the manner In which the load Is 
Introduced. The analysis presented herein assumes 
rigid body motion of the clamped edge; this assumption 
represents, a "worst case" condition. In any actual 
experiment there will almost certainly be some defor- 
mation and/or even slippage In the fixture. 

Quasi— Isotropic specimens respond uniaxially at 
locations at least 3/4 of the specimen width away from 
the edge; for cross-ply, the uniaxial range in 1.50 
width away from the edge. Since specimens tested In the 
IITRI fixture [1] have such short gage lengths, it may 
be concluded that a uniaxial .response can not be de- 
veloped in specimens* using this fixture. 


* An exception is unidirectional laminates with small 
values of . 


preceding page blank not riLMED 
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The constrained edge effect upon measured Young's 
modulus may be determined as follows. Let and E* 
denote, respectively, the actual modulus and experi- 
mentally determined modulus using strain gages at the 
location y along the center line. Thus 

E* - ^ — . (33) 

£3^(0, y) 

Combining (33) with Hooke's law (3) to eliminate the 
strain £3^(0 » 7)> we obtain 

a^ip,7) - V3^yOy(o,y) . ^ 

The measured strain e (o,y) will normally contain 
contributions from In-plane and out of plane bending. 

Since the stresses are odd function of y, the bending 
effects may be eliminated by using several gages and 
averaging' the results. 

Equation 34 has been evaluated for quasi-lsotroplc 
and cross-ply^ laminates at several locations y; the 
results are shown In Tables I6. and 17. Column (a) In each 
table indicates -the predicted experimental error If gages 
were placed -at y»0. .Similarly, column (b) in each table 


These results are base 
laminates, and v=0.900 




laminates . 




o 


shows the predicted error if gages were placed at 
Since, for each case, column (b) is closer to unity than 
column (a), placement of gages at y*±b/2 is a better lo- 
cation for strain gage placement.. In fact, calculations 
at other ▼alues Indicate y«±b/2 is the optimal location. 

Column (c) in Tables IS and 17 is based upon the 
assumption of’ three strain gages, two at the quarter 
points on one face and the third in the center of the 
opposite face. Clark and Lisagor [1] took extensive • 
measurements of graphite/epoxy using strain gages at pre- 
cisely^ these three points. Column (d) shows the experi- 
mental results based upon Clark and Lisagor 's original 
data^. It will be observed that comparison of the theo- 
retical results column Cc) with the experimental result 
column Cd) is exceptionally good for quasi-isotropic 
laminates. For cross-ply laminates. Table 17 shows a 

considerable discrepancy between predicted and actual 

# 

error^. The experimental results confirm the greater 
sensitivity of modulus to aspect ratio for the cross- 


1 Original stress-strain curve-i were available only for 
r*C.25r'0»50 and 1.0. For Q-I specimens, a "best-fit" 
straight line was constructed over the strain range 
E»0 to e *0.006. Average moduli for the three aspect 
ratios were 7.09x10 ksi, 6.71x10 kri, 6.39x10 ksi. 

The actual modulus was assumed to be 6.59x10 ksi for 
purposes of completing the column. 

2 Cclum Cd) cf Table 17 uses data as reported zy Clark and Llsaccr. 




























LIABLE 17* Predicted ezperlsental error of 

tauiii^'s modulus E for C-P laminate 
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plies, but not to the extent predicted. Presumably, the 
assumption of a rigid clamped edge Is not appropriate 
for short gage length, high Poisson's ratio- specimens 
In coo^reaslon using the UTRI fixture. 

A more plausible explanation* Is the rather hl^ 
stress levels near the constrained edge place the 
material well Into the non-linear ranve of behavior. 

[Note the high stress levels at bhe edges In Figs l6,17]. 
Consequently » It Is possible that the width of the 
specimen near the clamped edge expands non-linear ly, 
thereby greatly diminishing the constrained edge effect 

For completeness, we point out that Clark and 

Llsagor [l3 found that the modulus of unidirectional 
laminates was Independent of this Is consistent with 
the results of Chapter V. 

Although an explanation of compressive failure of 
composites was not one of the objectives of this study, 
some preliminary results are obtainable directly from 
the stress analysis. Failure theories for single plies 
may be applied directly to the stress distribution 
within each individual lamina. 


Delamination will occur when the interlaminar 


shear stresses t, and exceed the allowable loads 


for the epoxy. 


• c. tresses may be attrcximatelj 


Tnese shear 


obtained from the three-dimensional equilibriiua 
equations, i.e. 


1 

t^«C ?- 


aT** 


t^-C 


ax 

ar^ 


1 + AT. 


xz 


0 . 


ax 


ay 

aoj 4 


( 35 ) 


where the superscript i refers to the i-th lamina. At 
refers to the difference In value of shear stress across 
the 1-th .lamina, and t^ the. thickness of the i-th lamina. 
For small t^,' these shear stresses are very small, except 
where the in-plane stress exhibit large ^adients. 


Chapter VII 
COKCLUOINQ REMARKS 

Limitations of the Model 

Insofar as the problem is analyzed as generalized 

plane stress. It will not provide an exact solution to 

the three-dimensional elasticity problem* In particular, 

the third equation of equilibrium, (force-balance In the - 

z-dlrection) will not be satisfied D-51 . However, It is 

well known that the generalized plane stress solution is 

very close to the exact solution if the thickness of the 

• • * 

laminate la small compared to the other two dimensions. 

The linearity assumption Eqs. (3) Is a somewhat 
more serious limitation of this model. Compression 
tests of uniaxial C0]s high-strength graphite/epoxy 
laminates indicate linear behavior between load and 
axial compressive strain all the way' to fracture Ll3- 
Slnce the load is carried predominantly by the graphite 
fiber, it may be inferred that graphite responds linearly 
to compressive rupture. On the other hand, a cross-ply 
C±^5/?^5] s stacking of the same laminae produces a 
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non-linear behavior, particularly near failure. It Is 
Important to note that although the ultimate axial strain 
for eross«^l 7 . laminates exceeds the ultimate axial strain 
for unidirectional laminates by a factor of up to 3 C l3 « 
the maximum compressive fiber strain Is considerably 
lower for the cross-ply than for the uniaxial layup. 
Indeed, these cross-ply laminates fall due to delaml- 
. nation and not fracture [1]. It may be Inferred from 
the foregoing discussion, that the cross-ply laminate 
behaves non-Unearly because the epoxy exhibits non- 
linear behavior. Such non-Unear effects may also be 
observed from transverse strain measurements on 
unidirectional laminates. Ashton (16] reports varying 
values for Poisson's ratio during axial compression. tests 
on high-strength graphite/epoxy composites. The In- 
elastic behavior of composites was 'also Investigated ' 
by Poye- from the point view of- mlcromechanlcs £17]. 

The model Is very difficult to validate empiri- 
cally, since it is impossible to know the exact kinematic 
boundary conditions at the clamped edges. It is evident 
that an edge constrained to respond rigidly is the 
severest case that might be encountered. The results 
obtained in this study should therefore be viewed as 
the "worst possible case" . 


Treatment of the constrained edge effect due to 
out of plane bending, while of technical interest, is 
not studied in this work. Such effects are expected to 
be small in comparison to in-plane bendlns because the 
Poisson ratio is generally much smaller than Vxy wd 
the thickness of most laminates is very small compared 
to their width. Moreover, a study of these effects, 
would involve a considerably more complicated model. 

Thus the developed model should only be considered 
a first approximation to an accurate description. It 
may be used by the experimentalist to corroborate only 
the Inltal portion of the streca-straln compressive 
data. At the other end of the data curve it may be used 
only to suggest, rather than provide definitive 
explanations,, for different modes of failure. 
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